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SUMMARY 

A  mathematical  theory  of  the  effect  of  rapid  distortion  on  turbulence 
level,  originally  developed  by  Batchelor1  for  incompressible  flow,  has  been 
extended  to  account  for  compressibility.  Applied  to  a  wind-tunnel  contrac¬ 
tion  cone,  it  shows  that  the  effectiveness  of  the  contraction  in  reducing 
the  turbulence  level  and  non-uniformities  in  the  mean  flow  is  enhanced  as 
the  outlet  Mach  number  is  increased. 


Departmental  Reference:  Aero  3521 


Copyright 

© 

Controller  HMSO  London 
1982 


LIST  OF  CONTENTS 


1  INTRODUCTION 

2  EXTENSION  OF  THE  INCOMPRESSIBLE  THEORY 

3  ANALYSIS  FOR  A  GENERAL  DISTORTION 

4  APPLICATION  TO  A  WIND-TUNNEL  CONTRACTION  CC  'E 

5  CONCLUSIONS 

Appendix  A  Derivation  of  equation  (2-5) 

Appendix  B  The  form  of  the  energy  spectrum  tensor 

List  of  symbols 

References 

Illustrations 

Report  documentation  page 


Page 


3 

3 

5 


!0 
;  i 

12 

13 

1  /. 

2 

Figures  1-3 
inside  back  cover 


Accession  For 


NTIS  GRA&I 
DIIC  TAB 
Unannounced 
Justification- 


a 

□ 

□ 


By - 

Distribution/ 


Availability  Codes 


[Avail  and/or 
Special 


153 


3 


1  INTRODUCTION 

It  is  well  known  that  a  sudden  change  in  the  cross-sectional  area  of  a  stream  of 
fluid  has  a  substantial  effect  on  the  turbulent  kinetic  energy  as  well  as  on  the  dis¬ 
tribution  of  mean  velocity.  This  fact  is  used  in  the  design  of  wind  tunnels,  in  that  the 
working  section  is  generally  preceded  by  a  contraction  cone  in  order  to  obtain  a  near¬ 
uniform  flow  with  a  low  level  of  turbulence.  The  turbulence  intensity  is  not  in  fact 
reduced.  Rather,  it  increases  through  the  contraction,  but  the  acceleration  of  the  mean 
flow  is  correspondingly  greater  and  the  relative  turbulence  is  therefore  decreased. 

2 

A  mathematical  theory  of  such  rapid  distortions  has  been  developed.  Taylor 
showed  how  the  problem  may  be  linearised  and  looked  at  the  effect  of  a  contraction  on 
a  sinusoidal  disturbance,  taken  to  represent  one  component  of  the  turbulent  motion. 

1  3 

Batchelor  and  Batchelor  and  Proudman  then  examined  the  whole  spectrum  of  these  dis¬ 
turbances,  effectively  summing  over  all  contributions,  and  looked  at  arbitrary  as  well 
as  symmetric  contractions. 

However,  the  effect  of  compressibility  on  the  reduction  of  the  level  of  turbulence 
appears  not  to  have  been  looked  at  before  and  is  considered  in  this  Report.  It  emerges 
that  its  effect  is  beneficial.  As  the  Mach  number  of  the  flow  in  the  working  section 
increases,  the  contraction  ratio  needed  to  achieve  a  given  reduction  in  the  relative 
turbulence,  decreases.  It  will  be  shown  that  the  major  effect  of  compressibility  is  on 
the  acceleration  of  the  mean  flow.  In  contrast,  there  is  only  a  slight  increase  in  the 
turbulent  kinetic  energy. 

2  EXTENSION  OF  THE  INCOMPRESSIBLE  THEORY 

For  an  incompressible  flow,  Batchelor'  calculates  the  effect  of  rapid  distortion 
on  turbulent  kinetic  energy  by  linearising  the  equations  of  motion  and  relating  the 
vorticities  before  and  after  the  contraction.  The  same  procedure  may  be  repeated  for 
flows  in  which  compressibility  must  be  accounted  for. 

Consider  then  a  turbulent,  compressible  fluid,  in  steady,  uniform  mean  motion, 
suddenly  deformed.  Define  p  and  u  to  be  the  mean  density  and  velocity,  x  to  be 
the  position  of  a  general  element  of  fluid  and  p’  and  u'  to  be  the  turbulent  fluctu¬ 
ations  in  density  and  velocity.  The  implication  of  the  assumption  that  |u' j/(u|  is 
small,  which  allows  linearisation  of  the  problem  is  broadly  this:  the  distortion  is  so 
rapid  that  it  alone  affects  the  physical  attributes  of  a  particle  of  fluid,  such  as 
velocity,  density  and  ics  relative  displacement  from  any  other  particle. 

More  rigorously,  it  is  assumed  that  the  decay  of  turbulent  energy  occurs  on  a  much 

greater  time  scale  than  that  of  the  distortion  itself.  Thus,  the  convective  terms  in 

the  turbulent  kinetic  energy  equation  will  be  dominant.  Now,  the  rate  of  decay  of 

3  1 

turbulent  energy  is  proportional  to  |u' |  /£.  ,  (Batchelor  ) ,  where  l  is  the  length 
scale  of  the  turbulent  fluctuations  and  the  time  scale  for  decay  is  therefore  Jc/ju' [  , 
while  the  time  scale  of  the  distortion  is  of  order  D/|u[,  where  D  is  the  length  of 
the  distortion.  The  condition  under  which  rapid  distortion  theory  may  be  applied  is 
therefore 


(2-1) 


l 

D 


In  general,  l  is  small,  so  that  the  distortion  must  indeed  be  rapid,  if  the  linearisa¬ 
tion  shown  below  is  to  be  valid. 

When  (2-i)  holds,  there  is  insufficient  time  for  viscous  dissipation  to  affect  the 

turbulent  velocitv  and  Cauchv's  equation  can  be  used  to  link  the  vorticity  fluctuations 

4  . 

before  and  after  distortion.  This  equation,  given  by  Lamb  ,  is 


with 


and 


U)"(x)  =  oS„w'(a) 


12-2) 


Ox 


(2-3) 


3x. 

_ i 

3a. 

J 


(2-4) 


where  the  mean  density  of  the  fluid  is  px  after  distortion  and  ca  before  distortion, 
while  x  and  u"  are  the  position  and  vorticity  of  an  element  of  fluid  that  was 
originally  at  a  ,  with  vorticity  u'  .  This  analysis  does  not  consider  the  way  in 
which  vorticity  changes  during  the  contraction  and  assumes  uniform  mean  flow  before  and 
after  distortion. 

At  the  same  time,  because  the  Mach  number  of  the  turbulent  fluctuations  is  small, 
the  turbulent  part  of  the  flow  is  effectively  incompressible  and  conservation  of  mass 
implies  that 

p  V  .  u'  =  0  =  o  v  .  u"  .  (2-5) 

a-  -  x-  - 


The  derivation  of  this  equation  is  shown  in  Appendix  A. 

Also,  to  first  order  in  ju'|  ,  the  distortion  tensor,  is  related  only  to 

the  mean  velocity  and  is  therefore  uniform.  If  the  principle  axes  of  strain  are  used 
becomes  a  diagonal  matrix.  Analysis  of  the  problem  is  therefore  simplified 
and  the  continuity  equation  for  the  mean  flow  becomes,  in  Lagrangian  form 


e, e0e„ 
i  2  o 


a 


(2-6) 


where  e.  =  3..  . 

l  n 

Of  course,  the  arguments  presented  above  give  a  highly  idealised  description  of 
the  type  of  distortion  that  would  be  of  interest  in  practice,  since  other  effects  will 
generally  be  present.  For  instance,  if  the  distortion  is  caused  by  a  wind-tunnel  con¬ 
traction  cone  ,  there  will  be  enough  time  during  the  contraction  for  energy  exchange 
between  components  of  turbulence.  The  decay  of  turbulence  must  also  be  expected  to  be 
significant,  while  boundary  layers  forming  on  the  tunnel  walls  affect  the  mean  velocity 
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distribution  and  so  the  distortion  tensor  will  not  be  strictly  uniform.  The  work  which 
follows  should  therefore  be  thought  of  as  being  strictly  valid  only  in  the  limit  as  such 
secondary  effects  vanish.  It  is  nevertheless  useful,  because  it  gives  a  good  indication 
of  the  qualitative  effect  of  the  contraction  on  the  level  of  turbulence. 

3  ANALYSIS  FOR  A  GENERAL  DISTORTION 

Equations  (2-2)  to  (2-6)  form  a  basis  for  the  following  analysis.  For  the  sake  of 
generality,  arbitrary  values  of  ej,  e,  and  e^  are  assumed  at  first.  Specific  cases 
of  interest  will  be  considered  in  section  A. 

In  order  to  represent  the  properties  of  turbulence,  which  is  a  random  phenomenon, 
Batchelor's  approach  is  followed  and  Fourier  resolutions  of  the  velocity  fluctuations 
u*  and  u"  are  made,  thus: 

u'(a)  =  J  ei-‘-dzI  (k)  ,  u"(x)  =  Je1*"- dz"(x)  ,  (3-1) 


where  integration  is  carried  out  over  all  wave  number  space.  The  functions  dz 
(indices  '  and  "  will  be  dropped  unless  it  is  necessary  to  distinguish  between 
quantities  before  and  after  distortion)  are  statistically  random  and  provide  a  measure 
of  the  strength  of  the  turbulent  fluctuation  at  a  particular  wave  number.  Because  z(k) 
is  not  in  general  bounded,  the  integral  is  of  the  Fourier-Stieltjes  type.  Further  dis¬ 
cussion  of  the  properties  of  z  can  be  found  in  Batchelor1 . 


The  wave  number  vectors  k 
so  that 


and  x  are  related  via  the  distortion  tensor  S. . 
a  11 


i  =  1,2,3,  (3-2) 


while,  from  (2-6)  the  differential  volumes  dK]}  di^,  dtc^  and  dXj,  dx2>  dx3  are 
related  by  the  equation 


crdfCjdi^d^  =  dx,dx2dx3 


(3-3) 


A  relationship  between  dz'  and  dz"  may  also  be  found  on  substitution  of  (3-1)  into 
the  continuity  equation  (2-5)  which  yields 


k  ,  dz'  =  0  =  x  •  dz" 


(3-4) 


The  energy  spectrum  tensor  $,.(*) 


is  defined  by 


to 


u . (x)u . (x  +  r) 
i  -  j  -  - 


I 


> .  j  (x)  e1- ‘ -d*:  j  d<  2du  ^ 


(3-5) 


where 


/ 


<1*1^3 


< 


is  taken  over  all  wave  space. 


Batchelor 


1 


shows  that  this  definition 
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is  equivalent  to 


$.  •(<)  = 
ij  " 


lim 

dKjdK^dK^-K)' 


(..  Jzid2i  \ 

\°KldK2d«3/ 


(3-6) 


where  dz*  is  the  complex  conjugate  of  dz.  . 

Now  equation  (2-5)  implies  that,  with  the  summation  convention  used  for  repeated  indices. 


5u.(x  +  r) 

- 1_  =  o 


3r 


(3-7) 


and,  with  (3-5)  this  gives  that 


K  .  .  =  0 

J  IJ 


(3-8) 


or,  on  use  of  the  equivalent  expression  to  (3-5)  for  u^(x  -  r)u^(x)  , 


i  ij 


(3-9) 


As  a  consequence  of  the  definition,  it  can  be  shown  (see  Appendix  B)  that  the 

energy  spectrum  tensor  $£•  is  conpletely  defined  by  its  diagonal  elements.  Thus,  to 

f  3 

describe  ,  only  need  be  given  and  those  only  in  terms  of  ^ 

is  assumed  that  these  are  already  known.  To  find  such  a  relation,  take  the  curl  of 
equation  (2-2) ,  which  gives 

-vV(x)  =  -  (3-10) 


e.  32u’ (a) 


r  - 


"rsi'ipq  e  3a  3a 
n  s  p  s 


and,  with  (3-1),  implies  that 


2  e- 
X  dz"(x)  =  -  ae  .c.  —  k  k  dz'(tO 

r  -  rsi  lpq  eg  p  s  q  - 


(3-11) 


where  e.^  is  the  alternating  tensor. 

Expanded,  this  gives  three  equations  of  the  form 


i  /  e9  9  e9  9\  ®9  ) 

x)  =  -ajdz;(K)^-_K---,i-J  +  -1ci<2dz*  +^KiK3dz*J 


X~dz'J  (x) 

or,  using  equations  (2-6)  and  (3-2) 


(3-12) 


dz! (k) 


dzV(x)  = 

I  -  e 


<1  A 

“2  L 


x.dz! (tc) 
i  i  - 


X  el  “l  ei 


(3-13) 
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It  is  worth  noting  here  that  the  explicit  appearance  of  the  density  ratio  a  has  been 
eliminated  by  use  of  the  continuity  equation  (2-6) .  The  equation  is  then  identical  to 
that  obtained  by  Batchelor*  for  incompressible  flow.  However,  the  dependence  upon  a 
is  implicit  in  the  values  of  e^  . 

It  follows  from  equations  (3-3),  (3-6)  and  (3-17)  that  the  energy  spectrum  tensor 
is  given  by 


»re<*> 


o 

4 

X 


e . 

£  .£  .  — 1 

rsj  jpq  e 


e„  e  ,  —  kkkk  (<) 
xrav  vab  e  p  s  a  m  qb  - 
m  r 


(3-14) 


which,  after  some  manipulation,  and  use  of  the  relations  given  in  Appendix  B,  yields 
three  equations,  for  the  diagonal  elements  of  ,  of  the  form 


(3-15) 


It  is  now  possible  to  look  at  the  ratios  of  the  longitudinal  and  lateral  components 
of  the  turbulence  intensity,  given  by  p  and  v  , 


where 


-”2 

U1 

*2 

u. 


u  = 


_  X 


/*; 


jdX^XjdXj 


J*  ildKl 


dK2dK3 


-"2  ^  2 
u2  +u3 

2  *2 

u2  +  V 


=  v 


_  X_ 


I  (*22  +  *33>dXidX2dX3 


/  (*22  +  -?33)dKldK2dK3 


where  x^  is  chosen  as  the  longitudinal  coordinate, 
4  APPLICATION  TO  A  WIND-TUNNEL  CONTRACTION  CONE 


(3-16) 


(3-17) 


Consider  the  problem  of  a  symmetric  wind-tunnel  contraction  of  ratio  C  ,  as  in 
Fig  1 ,  and  say  a  cylindrical  element  of  fluid  has  radius  r  and  length  L  .  The  ratio 
of  the  cross-sectional  area  of  the  cylinder  to  that  of  the  tunnel  will  not  change  through 
the  contraction  cone,  so  that  the  radius  of  the  element  after  distortion  will  be  r/*^C 
and,  since  mass  is  conserved,  its  length  will  become  LC/o  .  From  this  argument,  it 
follows  that  the  values  of  e.  are  given  by 


A  simple  way  to  model  the  turbulence  upstream  of  the  contraction  is  to  assume  that 
it  is  isotropic  and  write,  following  Batchelor*, 


Substitution  of  these  values  and  the  values  of  e^  into  equation  (3-15)  and  the 
similar  expressions  for  ^  an<*  ^'33  laads,  after  further  manipulation  to  the  equations 

oE(k)(ic^  +  k“) 

ill  ^  \—  j/  /  / 


A"  ,  $ 

22  33 


f  2  2  A 

1  E(k)  Lv4  ^  K1K  0 


rx  +  — : _ 

.  2  4  _3 

4ittc  x  C 


which,  in  equations  (3-16)  and  (3-17),  yield 


-  [  ( sin^jA 

J  \fce); 

_0 _ 

J  sin^e  d9 


f  C  sin  e(  1  +  q>os2^ 
5  \  0^(9)  / 


sin  9(1  +  cos“9)  d9 


where  (k“  +  = 


x  sin  9  ,  *c  =  k  cos 


17 

"  j  [  dKjdK^dt*^  =  j  [  ^lrK~s^n  ®  dxd6 


/  2  \2 
f(9)  “  cos"9  +  sin“6  I 


. III! . Ill . I . . . . Ill . . . . 
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If  a  new  parameter,  a  ,  is  defined  by 


equations  (4-5)  and  (4-6)  reduce  to 


3c2  1  +  a2  ,  / 1  +  a\  J 

=  —  - —  *°s  e  —  -  — 

4C  2a  eV  '  a 


it  o 

3C  .  o”  1  1  -  a", 

-  '  *  J  72  -  TT-10 

C  ,2a  4a 


^i] 


the  total  turbulence  intensity  is  changed,  through  the  contraction,  by  a  factor 
which, because  the  turbulence  is  initially  isotropic,  is  given  by 


/  « - 
i.  I—  * 

u* .  /u.  - 
11/  11 


i(u  +  2v)  =  §  +  ~T-  lo8e(f^f) 


(4-10) 


The  effect  of  compressibility  through  the  distortion  can  now  be  seen-  When  a  -  1 

and  the  flow  is  incompressible,  the  expressions  (4-8)  to  (4-10)  reduce,  as  expected,  to 

those  given  by  3atchelor5.  Explicitly,  the  dependence  of  u  +  2v  on  o  (<  1  if  H  4  0) 

is  quadratic,  but  there  is  also  an  implicit  dependence  on  o  via  a  .  Fig  2  shows  a, 

u  and  v  plotted  against  a  range  of  values  of  the  contraction  ratio,  taking  the  outlet 

Mach  number  to  be  unity.  Clearly,  the  density  ratio  drops  quickly  and,  for  C  >  3  ,  is 

close  to  the  ratio  of  outlet  density  to  stagnation  density,  which  is  0-634.  The  value 

?  3 

of  a  ,  which  is  zero  at  C  =  1,  then  tends  to  unity  as  the  terra  o~/C  becomes  small. 
Correspondingly,  u  becomes  small,  showing  the  strong  effect  of  the  contraction  cone 
on  the  longitudinal  turbulence  intensity,  while  v  approximates  to  3 C/4  which  shows 
that  the  lateral  intensity  increases  through  the  contraction.  In  equation  (4-10)  it  can 
be  seen  that  for  C  >  3  the  lateral  dominates  the  longitudinal  contribution  and  that  the 
total  turbulence  intensity  increases  by  a  factor  which  is  approximately  C/2  ,  independ¬ 
ent  of  the  effect  of  compressibility. 

However,  the  factor  of  real  interest  is  that  which  applies  to  the  reduction  of 
percentage  turbulence,  ie  the  ratio  of  the  square  root  of  the  turbulence  intensity,  a 
measure  of  the  turbulent  velocity  fluctuations,  to  the  mean  velocity.  Because  mass  flux 
is  conserved  through  the  contraction  cone,  the  ratio  of  the  mean  velocities  before  and 
after  distortion  is  C/o  .  The  percentage  level  of  turbulence  is  therefore  reduced  by 
a  factor  F 


where 


r  ■ 


(4-11) 


and  compressibility  again  has  a  significant  effect.  For  large  values  or  C  ,  a  good 
approximation  to  F  is 


(4-12) 


that  is,  F  is  a  times  the  reduction  obtained  for  incompressible  flow.  This  shows 
that  the  contraction  ratio,  needed  to  obtain  a  particular  reduction  in  percentage 
turbulence,  decreases  as  the  outlet  Mach  number  increases.  Fig  3  compares  values  of  C 
for  various  percentage  turbulence  reductions  at  M  =  0,  M  =  0.5  and  M  =  1.  Further  results 
are  shown  in  Winter  and  Boyd"’ . 

Since  the  purpose  of  a  wind-tunnel  contraction  is  also  to  reduce  non-uniformities 
in  the  mean  flow,  it  is  valuable  to  look  at  the  effect  of  compressibility  in  this  con¬ 
text.  The  wave  number  vector  is  new  of  the  form  (0,  <^)  and,  if  the  contraction  is 

still  symmetric,  equation  (3-13)  gives,  with  equation  (3-4), 


dz’j  = 


dz?  =  C*dz! 
i  i 


(4-13) 


i  =  2,3 


which,  with  equations  (3-6),  (3-16)  and  (3-17)  show  that  the  percentage  level  of  non- 

2 ,  ? 

uniformity  in  the  streaawise  direction  is  reduced  by  a  factor  of  o  /C  while  the  angle 
of  flow  to  the  tunnel  centre  line,  caused  by  the  non-uniformities,  is  reduced  by  a  factor 
of  a/C*  .  In  both  cases,  compressibility  improves  the  effect  of  the  contraction. 

5  CONCLUSIONS 

(1)  The  theory  of  rapid  distortion,  derived  by  Batchelor1,  has  been  extended  to  account 
for  the  effect  of  compressibility. 

(2)  The  theory  has  been  applied  to  the  flow  through  a  contraction  and  it  emerges  that 
compressibility  has  only  a  small  effect  on  the  change  of  turbulence  intensity,  but  has 
a  significant  effect  on  the  reduction  of  percentage  total  turbulence. 

(3)  The  effect  of  compressibility  is  favourable.  The  contraction  ratio  needed  to 
achieve  a  given  reduction  in  the  level  of  turbulence  or  in  the  level  of  non-uniformity, 
decreases  as  the  outlet  Mach  number  increases. 


pPlilliPIIIPiPPPIIIIIII . IP . . . 
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Appendix  A 

DERIVATION  OF  EQUATION'  (2-5) 

It  is  assumed  that  the  flow  before  or  after  distortion  has  uniform  mean  density  o 
and  uniform  mean  velocity  U  in  the  longitudinal  direction  Xj  .  It  is  further  assumed 
that  the  turbulent  fluctuations  in  this  region  are,  to  first  order  in  0([u' |/U) ,  simply 
convected  by  the  fluid,  so  that 


du' 

dt“ 


0 


dp' 

dt 


(A-l) 


In  fact,  this  assumes  that  there  is  no  mechanism  present,  to  first  order,  by  which 
turbulence  is  created  or  destroyed  and  implies  that  the  fluctuation  in  pressure,  p  , 
is  0(iu' |~/D“)  - 

Now,  the  principle  of  continuity  of  mass  gives 


')  + 


(o  +  p'Hu 


=  0 


(A-2) 


or,  to  first  order  in  ]u'|/U  ,  since  P  and  U  are  constant. 


3P* 

3t 


0  - 


(A-3) 


Combination  of  this  equation  with  equation  (A-l)  gives  the  required  result,  equation  (2-5) 
or 


3uj  3U; 

o  - +  — —  + 

rxi  3x2 


0 


(A-4) 


which  states  that  the  turbulent  velocity  fluctuations  are  incompressible  to  first  order* 
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Appendix  B 

THE  FORM  OF  THE  ENERGY  SPECTRUM  TENSOR 

The  energy  spectrum  tensor  ,  as  defined  in  equation  (3-6),  is  clearly  hermitian, 

so  its  diagonal  elements  are  all  real,  while  its  off-diagonal  elements  are  related  by 
three  expressions  like 

$.  .  =  <j>*.  j  t  i  •  ( B—  1 ) 

ij  Ji 


There  are,  therefore,  only  nine  independent,  non-zero  quantities  in  the  set 


{re$.  .  ,im4>.  . } . 
1J  ij  i 


Equations  (3-8)  and  (3-9)  together  with  relations,  which  follow 


i.J 


from  equation  (3-6),  of  the  form 


=  4>.  - 4> -  -  ;  i  ^  i  ,  both  fixed 

ii  JJ  i]  ]i 

.  (B-2) 

apparently  give  nine  equations  relating  the  members  of  this  set.  However,  on  closer 
inspection,  it  emerges  that  only  six  of  these  nine  are  independent  and  these  are:  from 
equations  (3-8)  and  (3-9),  with  (B—  1 ) 


K|te<J> 

iCjre^ 

«,re4> 


11 

12 

13 


+ 


+ 


+ 


K„re$ 

£. 

«2re?> 

K2re<ji 

<2^ 
k  j  im$ 


12 

22 

23 

12 

12 


+ 


+ 


+ 


+ 


+ 


x^reifi 

x^reij) 

x3reiji 

x^imij) 

K^imiji 


13 

23 

33 

13 

23 


and  from  equation  (B-2)  with  (B— 1 ) 


'11*22 


=  (re^12)2  + 


(im<J>  j  2) 


(B-3) 


(B-4) 


Three  further  equations  only  are  needed  now  to  define  ^  completely  and  are 
provided  by  equations  giving  the  diagonal  elements,  ^j,  <J>22  and  $33  .  Thus,  the 
energy  spectrum  -ensor  <J> . .  is  completely  determined  by  its  diagonal  elements. 


153 
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LIST  OF  SYMBOLS 


position  before  distortion  -  (aj,a2,a3) 


u',  u" 


contraction  ratio 

length  scale  of  distortion 

energy  function  (equation  (4-2)) 

factor  of  reduction  of  percentage  turbulence  through  distortion 

/-T 

turbulent  length  scale 

length  r-f  cylindrical  element  of  fluid  before  distortion 

outlet  Mach  number 

pressure 

radius  of  cylindrical  element  of  fluid  before  distortion 

position  vector  =  (r|>r2,r3^ 

distortion  tensor 

mean  velocity  vector  =  (U^U^U^) 

turbulent  velocity  fluctuations 

position  after  distortion  =  (x^x^.x^) 

measure  of  the  strength  'f  turbulent  velocity  fluctuations  =  (dZj ,dz2 ,d23) 
complex  conjugate  of  dz  _ _ _ 

r  2  3 

defined  in  equation  (4-7)  =  v  1  -  (o  /C  ) 

the  alternating  tensor  =0  i=j»j=k  or  k=i 

=  1  ijk  are  cyclic 

=  -1  ijk  are  anti-cyclic 

wave  number  vector  before  distortion  =  ^],|J2,K3^ 

ratio  of  longitudinal  components  of  turbulence  intensities  before  and  after 
distortion 

ratio  of  lateral  components  of  turbulence  intensities  before  and  after  distortion 
mean  density 

turbulent  fluctuation  of  density 

ratio  of  mean  densities  before  and  after  contraction 
energy  spectrum  tensor 

wave  number  vector  after  distortion  =  (Xj>X2*X2^ 
vorticity 
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